Exact 1-D Model for Coherent Synchrotron Radiation with Shielding and Bunch 
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Coherent Synchrotron Radiation has been studied effectively using a f-dimensional model for 
the charge distribution in the realm of small angle approximations and high energies. Here we 
use Jefimenko's form of Maxwell's equations, without such approximations, to calculate the exact 
wake-fields due to this effect in multiple bends and drifts. It has been shown before that the 
influence of a drift can propagate well into a subsequent bend. We show, for reasonable parameters, 
that the influence of a previous bend can also propagate well into a subsequent bend, and that 
this is especially important at the beginning of a bend. Shielding by conducting parallel plates is 
simulated using the image charge method. We extend the formalism to situations with compressing 
and decompressing distributions, and conclude that simpler approximations to bunch compression 
usually overestimates the effect. Additionally, an exact formula for the coherent power radiated by 
a Gaussian bunch is derived by considering the coherent synchrotron radiation spectrum, and is 
used to check the accuracy of wake-field calculations. 



I. INTRODUCTION 

Coherent Synchrotron Radiation (CSR) is an impor- 
tant detrimental effect in modern particle accelerators 
with high bunch charges and short bunch lengths. It 
is a collective phenomenon where the energy radiated 
at wavelengths on the order of the bunch length is en- 
hanced by the number of charges in the bunch. This CSR 
field can subsequently affect particle motion. A compre- 
hensive collection of papers and references regarding this 
subject is found in Ref. 

CSR is difficult to model using discrete particles ex- 
actly because the problem scales with the number of 
particles N as A^ 2 . To simplify this, the 1-dimensional 
model projects the transverse particle density onto the 
longitudinal dimension. Formulas for the CSR field from 
this line charge are then used to calculate forces on each 
particle and then propagate the full bunch distribution. 
While this makes the calculation tractable, the electro- 
magnetic fields on the world-sheet of this charged line 
are singular. Pioneering efforts described in Refs. @, Q 
circumvent this problem by examining the non-singular 
terms only. In more detail, Ref. 0] regularizes the lon- 
gitudinal force between two charges traveling on the arc 
of a circle by subtracting off the Coulomb force, calcu- 
lated as if the charges traveled on a straight line, from 
the force calculated using Lienard-Wiechert fields. The 
result is an always finite CSR force. 

This paper uses the less widely known Jefimcnko forms 
of Maxwell's equations @, which allow one to calcu- 
late electromagnetic fields by directly using the evolving 
charge and current densities, and which internally incor- 
porate all retardation effects. These equations are re- 
lated to forms used in Refs. 0, This is in contrast 
to the usual Lienard-Wiechert approach, which gives 
fields due to charges at their retarded times t' and po- 
sitions x(t'), and one must invert equations of the form 
t — t'= |x(t') — x Q |/c for the retarded time t', where x D is 
an observation point at a later time t and c is the speed 



of light. While this latter method has proven useful in 
deriving equations for (incoherent) synchrotron radiation 
of single particles, the former is found to be useful for the 
coherent fields of particle distributions. 



II. EXACT ID APPROACHES TO STEADY 
STATE CSR 

In general, for given charge and current densities 
p(x,t) and J(x, f) at position x and time t, the electric 
field E(x, t) can be calculated using Jefimenko's form of 
Maxwell's equations [5[ 
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t'=t-r/c 



in which r = x — x', r = |jr||, eo is the vacuum permit- 
tivity and t' is the retarded time. In this formulation, 
the retarded points x' and times t' are independent vari- 
ables, so there are no functions that need to be inverted. 
Therefore, if one knows p, p, and J at all points in space 
x' and times t' < i, with a dot denoting the time deriva- 
tive, then this formula gives the electric field by direct 
integration. 

Now consider a line charge distribution, which follows 
a path X(s) parameterized by distance s, has a unit tan- 
gent u(s) = dX.(s)/ds, and moves with constant speed 
/3c along this path. A bunch with total charge Q and 
normalized line density A therefore has one-dimensional 
charge density and current 



p(s, t) — Q X(s — s;, — j3 ct), 
J(s, t) = Q/3cu(s) A(s — si, - 



pet), 



(2) 



where Sb is the location of the bunch center at time t = 0. 

The rate of energy change per unit length of an ele- 
mentary charge q at position s is d£ /ds — q u(s) • E(s, t). 
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Functions of this type are called wake-fields. Using 
Eq. |T]) with the one-dimensional bunch in Eq. (0) gives 



— (s,t) = Nr c mc 2 
as 



u(s) -r(s,s') 
r(s, s') 3 



A(s r ) 



(3) 
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u( S ) 



r(s, s') 2 
with the definitions 



A'(s r ) 



a2 u(s)-u(s') 

P 7 7^ — 

r(s, s'j 



A'(s r ) 



s r = s' - s +fir{s, s), 

s = S(, + /3 ci, 
r(s,s') = X(s)-X(s'), 
r(s,s')= ||r(s,s')||, 



(4) 
(5) 
(6) 
(7) 



where iV = Q/q is the number of elementary particles 
with mass m and classical radius r c = q 2 / (47reomc 2 ) , 
and the prime on A indicates a derivative of this func- 
tion with respect to its argument, i.e. A'(x) = dX/dx. 
Additionally, sq is the center of the bunch at time t, and 
this is the only place where the time dependence appears. 
The integrand is thus the contribution to the wake-field 
due to particles between the retarded positions s' and 
s' + ds', with N A(s r ) being the charge density at re- 
tarded position s' and retarded time t' . 

Unfortunately the integral in Eq. Q diverges as s - 
s' — > 0, which is a consequence of the one-dimensional 
line charge model. This problem can be alleviated by 
using the regularization procedure originating in Saldin 
et al. [H , where the electric field E is split into two parts 

E = E CSR + E sc . (8) 

The space charge (SC) part is the electric field of a line 
charge moving on a straight path, 



EscM) 



__Q_ 
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u(s) 



rA(s ; ) 



(9) 



\s — s 

with si = J — so + (3 \ s 
parts, simplifying to 



rA'( S/ )+/? 2 



s\, which can be integrated by 



Esc(M) 



Qu(s) 
47reo7 2 



ds 



A'(s - s + (3\s - s|) 



(10) 



It will turn out to be useful to change variables in this 
expression, so that when combined with Eq. the func- 
tion A' can be factored. This can be done by setting 
J — so + /3 \s — sj = s' — sq + 13 r(s, s'), with the con- 
vention that sgn(s — s) = sgn(s — s'). Noting that 
dr(s, s')/ds' = — r(s, s') ■ u(s')/r(s, s'), this leads to 

T , , . 1 + (3 sgn(s' - s) 



= sgn(s' — s) 



ds = 



1 



s — s' — (3r(s, s') ' 
f3r(s,s')-u(s')/r(s,s') 



1 + (3 sgn(s' — s) 



ds', 



(11) 
(12) 



so that 

Eso(M) 



u(s) 
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ds' sgn(s' — s)A'(s r ) 



(13) 



l-/3r(s,s') ■ u(s')/r(s,s') 
s — s' — /3r(s, s') 



The resulting wake-field due to E CSR , called the CSR- 
wake, is 



d£cSR 

ds 



= qu(s)-[E(s,t)-E sc (s,t)} 



(14) 



This expression is finite, and shown in Ref. [J] to correctly 
account for the coherent energy loss due to synchrotron 
radiation. 

The approach here is to be contrasted with the conven- 
tional one taken in the literature using Lienard-Wiechert 
formulas. In terms of the quantities above, the electric 
field at position s due to a charge q at retarded time 
t' = t — r(s, s')/c and retarded position s' is 



E L „(s,s') = 



1 



Pru(s') 



4tt€ [ 7 2 [r-/?r • u(s')] 

r x {[r- j3ru(s')] x ^ 2 u'(s')} 
[r-/3r -u(s')] 3 



(15) 



with r as in Eq. ^ suppressing the arguments. There- 
fore, the electric field at s due to a charge p(st, t) dst be- 
tween St and St + dst, as in Eq. is found by inverting 
St = s' + /3r(s, s') for s' and using Eq. (fT5")) . This is often 
impossible to do analytically, but fortunately for a distri- 
bution of charges the inversion can be circumvented by 
changing variables. Because dr(s,s')/ds' = — r • u(s')/r 
from before, the charge is 



p(s t ,t) ds t —Q A(s' — Sf, — (3 ct + (3r) 



r • u(s 
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r 

and the total electric field is 

oo 

E(s,*)= J ds t E LW (s,s'(s t )) p(s t ,t) 

— oo 
oo 

r • u(s') 



ds', 



(16) 



(17) 



Q / ds' 



1-/3- 



E LW (s,s')A(s r ). (18) 



We can use this to verify that E sc computed this way 
agrees with the result using the Jefimenko approach. For 
the SC field one has r(s) = (s — s')u(s), u(s') = u(s), 
and u'(s) = 0, giving 

oo 

Esc(M) = ^ / ds ,/3 + sgn(s-s') 

y ' 47reo J (s-s') 2 (19) 

— oo 

x \ {s' -s +/3\s-s'\). 
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FIG. 1: Geometry for a single bend. The variable s parame- 
terizes the curve with radius 1/k. The coordinates are X(s), 
and the unit tangent vector is u(s). 



Equation ([19]) agrees with Eq. (|10[) when integrated by 
parts because, for so = 0, / (/3+sgn(s— s'))(s— s')~ 2 ds' = 

~-A(s'+/3(s-s')sgn(s- 



(/3+sgn(s-s'))(s-s') , and 

s')) = -(l-/3sgn(s-s'))A'(s' + /3|s-s'|), and similarly 
for all Sq. 



III. SINGLE BENDING MAGNET 

Now we apply Eq. to the geometry of an arc of 
a circle of curvature n and length B, shown in Fig. [T] 
Set s = at the entrance of the bend so that = ks 
is the angle into the bend. In terms of fixed Cartesian 
unit vectors e a and e&, the path coordinates and tangent 
vector are 

X(s) = sin (k s) e a — k^ 1 [1 — cos (k s)] e&, (20) 
u(s) = cos (k s) e a — sin (k s) e&. (21) 

Consider a bunch with its center at angle 9q = ksq, 
and a test particle at angle 9. The contribution to Eq. © 
of this finite arc is 



d£ 
ds 



(s) 



Nr r mc 2 



d9' 



sin a 

{nr a ) z 



n\(s a ) 



(22) 



with s 



„ sin a , , . . „o cos a w . . 
~h A ( s «) + /? A'(s a ) 

_1 and the following definitions: 
a = 0-0', 

s a = - (0 - 0o - a) + /9r a , 

K 
1 



K 



2 cos a. 



(23) 
(24) 

(25) 



Thus a is the angle between the test particle and the 
retarded source particle, and is positive when the former 



is ahead of the latter. The first term of Eq. (j22f can 
be integrated by parts because 0(2 — 2 cos a) -1 / 2 / 89' = 
sin (a) (2 — 2cosa)~ 3 / 2 , and the wake greatly simplifies 
to 



dS 
ds 



—Nr c mc 



-K X(s a ) 



\J2 — 2 cos a 



«=-(« B-b 



da 



/3 2 cos (a) — 1 , 



(26) 



-(kB-0) 



\J2 — 2 cos a 



In terms of the variable a, the space charge term in 
Eq. (fl3|) can be split as 



d£ s 



-(kB-6) 



ds 



Nr c mc 



J dal sc (a) 



(27) 



+ / da/ sc (a)+ / dal sc (a)}, 
-{k.b-6) e 

with the integrand 

j3 sin (a) 

i 

«cm n 

Isc(a) 



sgncv 1 y/2 - 2 cos I 



7 2 a — (3 \J2~2 cos a 



A'(Sa), (28) 



so that the contribution of the bend to the CSR-wake is 



d£ c 



ds 



(s) 



-kA(s ) 



= Nr c mc 2 \ —j^^=^^= 
B I v 2 — 2 cos a 



-(kB- 



dcv A'(s a ) 

-(kB-0) 



(i 2 cos (a) - 1 
V2 — 2 cos a 

/3sin(a) (29) 



+ 



sgn(a) y/2 - 2 cos a 
7 2 a — /3 \J2 - 2 cos a 



-(« b-0) 



y da 7 S0 (a) ~ J da I sc (a) I . 



A. Steady State 

In the practical environment of a particle accelerator 
with a bunched beam, one is typically only concerned 
with electric fields around the bunch center. Due to the 
rotational symmetry, there will be an angle into a bend- 
ing magnet beyond which the CSR-wake, relative to the 
bunch center, does not change. Note that in Eq. (|24[) the 
quantity z = k^ 1 (9 — 9q) is the distance along the path 
ahead of the bunch center, and define the extent of the 
bunch If, = z + — z_ , where z + is the head particle coor- 
dinate, and z_ is the tail particle coordinate. Henceforth 
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FIG. 2: The inverse of Eq. (|33[l for positive a at various 
energies. The dashed green curve (7 — > 00) is a = (24 re A) 1//3 , 
the inverse of Eq. (|35j) . 



the symbol 2 will refer to the longitudinal coordinate rel- 
ative to the bunch center: z = s — sq. The particle at 
z + is affected by a particle at z_ at retarded angle a max 
found by inverting 



nlh 



(3\/2 — 2 cos a u 



(30) 



Similarly, a particle at z_ is affected by a particle at z + 
at retarded angle a m i n found by inverting 



- Klh = 



- /V2 - 2 cost 



When the bunch center is at an angle #0 > Qn 



(31) 



only particles within the bend affect the wake-field. The 
"steady-state" (s.s.) CSR-wake is then 



Wcsr(z) = Nr c mc 2 



da 



(3 2 cos (a) — 1 
y/2 — 2 cos a 



_^ (3 sin (a) 
sgn(a) V2 - 2 cos a 



(32) 



a- 



/3 y/2 - 2 cos a 



A' (z-A(a)) 



where 



A(a) = K' 1 (a - /3 V2 - 2 cos a) (33) 
is the distance behind the test particle at z. The notation 



W CSR (z) 



'-(s + z) 



(34) 



is used to refer to the CSR-wake immediately surround- 
ing the bunch center at sq. 

In the ultra-relativistic approximation (/3 — ► 1) with a 
small normalized bunch length kI^ <C 1, and thus a«l, 
the steady-state formula in Eq. (|32[) greatly simplifies. 
The I/7 2 term in Eq. (j2"5j) puts J sc — > 0, and the term 
in the integrand (/3 2 cos (a) - l)/(2 |sin (a/2)|) « -|a|/2. 
The function A (a) for 7 — > 00 is approximately 



A(a) 



a 3 /(24k) for a > 



2c*/k 



for a < 



(35) 
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FIG. 3: The steady-state CSR-wake for various relativistic 7 
using Eq. (I32|l . compared to Eq. (|36p plotted as green. Here 
re a = 3 x 10~ 5 for a Gaussian bunch, represented in light 
blue. 
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FIG. 4: The steady-state CSR-wake for various relativistic 7 
due only to particles ahead of the test particle, i.e. negative a 
in Eq. 1|32[) . and the second term in the integrand of Eq. (|36[1 . 
A Gaussian bunch is used, with recr = 3 x 10 -5 , and the 
wake has been scaled by (recr)" 4 / 3 . Compared with Fig. [3] 
this demonstrates that the contribution to the CSR-wake of 
particles ahead of the test particle is insignificant compared 
to those behind. 



Figure [2] plots the inverse of Eq. (f3"3"]l for positive a 
and various energies. One sees that the approximation 
in Eq. (|35[) is increasingly good for higher energies, but 
greatly overestimates a at the smallest distances. Chang- 
ing variables using Eq. ([35|) , the ultra-relativistic steady- 
state CSR-wake is 



Wcsr (z) 

7 — >oo 







2 \'{z - A) 

(3kA)V3 



K A ... . . 
+ -5- A (2 + A) 



(36) 



The first term in this integral is derived by an alternate 
method in Ref. [4] . The scaling here is apparent by writ- 
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FIG. 5: The ratio of the average energy of a Gaussian bunch 
using the exact Eq. (|32p to that using approximate Eq. (|36p 
in a practical range of the parameters 7 and kg. 
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FIG. 6: Similar to Fig. [5] but with the ratio of the average 
energy of Eq. (132[l to Eq. (|4ip. showing that the latter is an 
excellent approximation at relatively low energies. 



ing the distribution in the normalized form 



l~/z-A 

X(z - A) = -A 

a \ <j 



X'(z - A) = \\' ( : A 



(37) 
(38) 



where a 2 is the variance of A, so that A has unit variance. 
Also using normalized z = z/a and A = A/er gives 

( Kcr )2/3 

Wcsr (za) = -Nr c mc 2 - ^ — 



7 — >oo 
h/cr 

x/dA 

II 



a- 



^^ + (-) 4/3 ^ + A) 
(3A)V3 8 



1 (39) 



Now one can see that the particles in front of the test 
particle, represented in the last term in the integrand, 
influence the wake by roughly a factor of (kct) 4 / 3 less 
than particles behind, and that the primary contribution 
to the CSR-wake scales with the factor in front of the 
integral in Eq. (|39|) . However, it is interesting to note 
that even as 7 — > 00, where a charge radiates infinitely 
more power in the forward direction than the backward 
direction, there is still a finite CSR force from particles 
ahead of the test particle. In light of the primary scaling, 
we define a characteristic CSR energy change per unit 
length as 



W = Nr c mc 



(k a 



1 2/3 



(40) 



The ultra-relativistic approximation in Eq. (|36[) is com- 
pared to the exact formula Eq. ([3^)1 in Fig. [3] for var- 
ious energies and a particular value of n a. One sees 



that Eq. (|5rj)) represents the largest possible effect. The 
CSR-wake due only to particles in front of the test par- 
ticle is shown in Fig. IU emphasizing again that these 
forward particles contribute only a small amount to the 
total CSR-wake. 



Neglecting the contribution due to forward particles, 
the ultra-relativistic stead-state CSR-wake in Eq. ([39]) 
scales with Wo and depends only on the shape of A. Fac- 
toring out Wq from the exact steady-state CSR-wake in 
Eq. (|32l) . the exact result additionally depends on 7 and 
k 0. Therefore, to quantify the appropriateness of the 
ultra-relativistic approximation, the ratio of the average 
energy lost (per unit length) of a Gaussian bunch using 
the exact Eq. |32|) to that using approximate Eq. |36|) 
is shown in Fig. [5] for a practical range of these parame- 
ters. At a given energy, one sees that Eq. (|3"fjj) is a good 
approximation for the relatively long bunches. This can 
be understood from Fig. [21 because the approximation 
in Eq. ()35|) has a relative error for a finite energy that 
diverges for small a. 



A systematic method for calculating the CSR-wake us- 
ing Lienard-Wicchcrt formulas in the small angle, rela- 
tivistic approximations has been developed in Ref. 0] for 
arbitrary combinations of drifts and bends. Using the 
corresponding equation in Ref. [7j for the geometry of a 
bend, and the appropriate Jacobian factor, the steady- 
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FIG. 7: The free space and shielded CSR-wakes in a bend. 
The contributions to the shielded wake of individual image 
bunches are shown in red and blue. A lGeV Gaussian bunch 
with a — 0.3mm is used in a bend of radius = 10.0m. 
The shielding height H = 2cm. 




FIG. 8: The steady-state retarded distribution A re t,s.s. (z'; z±) 
for various test particles Zi in Eq. (|5ip using a Gaussian bunch 
with standard deviation a = 0.3mm and energy lGeV, in a 
magnet of bending radius xT 1 = 10.0m. 



state CSR-wake to second order in a and I/7 is 

1 



Wcsr (z) = —Nr c mc 

SHMS,s.s. v ' 



3 J d " 



2„,2 



x A' 



a 

■- 

2 + 7 2 a 2 1 

cv + 7 2 a 3 /4 a/2 + 7 2 a 3 /24 

a a 3 



(41) 



2 7 2 24 



Compared to Eq. (136p . this expression is a significantly 
better approximation of Eq. (|3"2")l for low 7 and a practical 
range of feu, shown in Fig. [6] 



B. Shielding by Parallel Plates 

The presence of a conducting beam chamber can have 
a strong effect on the CSR wake-field. For a rectangular 
cross section, it has been observed that the dominant 
effect comes from the smaller of the height and width 
(see, for example, Ref. Q). If particle trajectories are 
planar, then a finite chamber height can be represented 
by infinite parallel plates. In such a geometry, CSR wake- 
fields can be calculated using the image charge method. 

The kick due to a single image bunch at height h is 
easily adapted from Eq. ([3]) as 



d£ , , , „ T 9 f , , I u(s) ■ r , . , 
-(M,M = Nr c mJ J As' (? , 2 ^ 2)3/2 A( Sft ) 



(42) 



' 2 u(s) ■ u(s') _ u( 8 ) -r 

P ( r 2 + /j2)l/2 ^ r 2 + /j2 



A'(s fc ) 



with the argument 



Sh 



s b - (3ct + (3\/r 2 + h 2 , 



(43) 



and with r and u retaining their meaning from Eq. ([3]). 
Parallel plates require an image bunch for each plate, and 
an image bunch for each of those, ad infinitum. For the 
real bunch with orbit midway between plates separated 
by a distance H, symmetry gives the total image kick 



images 

ds 



d£ 



(s,t)= (-l) n ^(s,t,nH 



n— — oo 



2 E(- 1 )"x( s ' i ' ni/ )- 



ds 



(44) 



(45) 



If the real bunch has a vertical offset V, the total image 
kick is modified to 



d£ lT , 



ds 



E 



n— — oo.n^O 
even 



d£ 
ds 



s,t,nH) 



d£ 



(46) 



Ts (s,t,nH-2V). 



n— — CG 
odd 



In a bend, the contribution of the image bunches to the 
CSR-wake within the bend, following Eq. ([26]) . is 



d£h 



ds 



-(*) 



E 2 M 



Nr c mc 2 

n \ k A(s QjI j) 



-( K L m -0) 



(47) 



da 

-(«L m -e) 

with the definitions 



P 2 cos (a) — 1 , 



A (Sa,n) 



r Qj „ = y/2 — 2 cos a + (nnH) 2 , 
S a ,n = n^ 1 {0 - 6»o - a + j3 r a , n ) ■ 



(48) 
(49) 



Zt = 3cr, 2cr, lcr, 0, -ltr, -2o% -3cr 
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FIG. 9: The same as Fig. [HI along with image charges at a heights nH — n x 2cm (not to scale), which are approximately at 
heights n x 67a, and calculated using Eq. (|52[1 . 



Notice that the integrands do not need to be regularized 
by the SC term, because they are always finite due to the 
always positive factor (nnH) 2 . 

Due to the infinite number of image layers needed, 
a finite bend can never be exactly in the steady-state. 
However, due to their increased distances and angles, the 
relevant contribution image number n will be negligible 
beyond some maximum image number. This point is il- 
lustrated in Fig. [71 where the contributions to the CSR- 
wake of five individual images are shown along with their 
sum with the free space wake, to give the total shielded 
wake. 



C. Retarded Bunch Visualization 

For a given particle at time t within the bunch, it 
is evident that the retarded bunch density can be very 
distorted relative to the actual bunch density. From 
Eq. (fT6|) . the retarded bunch density at position s' as 
seen by a particle at position s is 



A rct (s';s) = X(s'-s b -/3ct+l3r) 1 



P V -^-) ■ (50) 



This retarded density is illustrated in Fig. [S] for a Gaus- 
sian bunch distribution for various test particles. There 
one sees that the density in front of the test particle is 
compressed to roughly cr/(l + /3) w a/2, concentrated in 
an apparent spike at the right of the plot. The density 
behind the test particle occupies the majority of the plot. 
While it may seem that the curves shown are Gaussian 
in form, this is only true for the left sides of the curves; 
the right sides have been extended and diluted due to 
the Jacobian factor in Eq. (|50| . Similarly, the retarded 
density of an image bunch at height h is 



A rot (s', h; s) 



1-/3 



r • u(s') 



Vr 2 + /i 2 



(52) 



x A(s' - s b - 0ct + l3\/r 2 + h 2 ). 



Figure [5] shows the retarded densities for a Gaussian 
bunch and several image bunches within a bend. In this 
example, the first and second image bunches as seen by 
particles in the rear of the bunch are actually closer than 
the real retarded bunch. 



IV. MULTIPLE BENDS AND DRIFTS 



In the steady-state, the geometry of a bend can be used 
in Eq. ([50)) . Moving to coordinates relative to the bunch 
center, the steady-state density seen by a test particle at 
Zt within the bunch as a function of z' is 



\ rc t (z ; Zt) 

s.s. 



1-/3- 



sin(/c(z t - z')) 



y/2 - 2cos(ft(z t - z')) 
A(V + /3 K - 1 v /2 _ 2 cos(k(z 4 - z')) ) . 



(51) 



A. CSR in Multiple Bends 

In this section the general formula Eq. ([3|) , regularized 
by Eq. (flT)]) . is applied to the geometry of multiply con- 
nected bends and drifts. Shielding by conducting parallel 
plates is added as in Eq. (|45)) . It has been seen in Eq. ([39]) 
that the primary contribution to the CSR- wake in a bend 
is due to particles behind the test particle, so for brevity 
the path is given behind the test particle only. 
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FIG. 10: Geometry for bends and drifts. The variable s pa- 
rameterizes the path X(s), with s = at the beginning of 
element Bi. The names Bi, D\, etc., also serve to indicate 
the element length. The dashed line is for a prior bend with 
negative curvature. 



Let the bunch center be at length so inside bend 1 of 
length B\ and positive curvature k\ , preceded by drift 1 
of length Di, preceded by bend 2 of length B2 and cur- 
vature K2 7^ 0, as shown in Fig.[l0j A drift follows bend 
1, referred to as Dq. A negative curvature K2 signifies a 
bend in the opposite direction of bend 1. With s = 
located at the beginning of bend 1 , the path coordinates 
are 

Xfl (s) for s > B\ 
X Bl (s) forO<s<Si 
s e a for —D\ < s < 
X Bj (s) far*<-Di 



X(s) 



(53) 



where the paths in the individual elements are 
sin (k\Bi) 



Xdo(s) 



+ (s — Bi) cos (k.\Bi) 



1 



(cos (kiBi) — 1) — (s — Bi) sin (kiBi) 



(54) 



e b , 



sin (ki s) [1 - cos (ki s)] 
X Bl (s) = e a e 6 , (55) 



X B2 (s) 



sin (k 2 Di + k 2 s) 



K2 



1 — cos (k 2 Di + K2 s) 



(56) 



K 2 



The tangent vector is then 

cos (mBi) e a 

— sin (k\Bi) eb 
cos (k% s) e a 
u(s) — ^ — sin [k\ s) Sb 
le a 



for s > Bi 

for < s < B 1 
for -Di < s < 



(57) 



cos (k 2 -Di + «2S) e a , / n 
for s < — Di 

— sin (k 2 Di + k 2 s) e~b 
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FIG. 11: The average energy loss and energy spread in- 
duced, per unit length, of the CSR-wake for a Gaussian bunch 
through the length of a bend in free space as well as between 
parallel plates with H = 2cm. Solid lines have Di — lm, 
while dashed lines have D\ — > 00. Parameters used are 
Kj -1 = k^ 1 = 10m, <t z = 0.3mm, with an energy of 1 GeV. 



Straightforward calculation gives the total CSR-wake at 
position s in the bend (0 < s < Bi) due to these different 
sections of the path 



d£ c 



ds 



(0 < s < 5i) = 



d£ a 



ds 



ds 



ds 



(58) 



+ 



with Bi, D±, and B2 signifying the contributions from 
bend 1, drift 1, and bend 2, respectively. Due to their 
length, these terms are written out in Appendix [A] 

A visualization of the retarded bunch and images of 
this geometry, similar to Fig. [HI is shown in Fig. [T2J Even 
though the bunch has progressed 50cm into bend 1, it sees 
much of the retarded bunch inside bend 2, especially for 
test particles z t in the front of the bunch. 

Two principal effects of CSR on the bunch distribution 
are a loss of energy and an increase in energy spread. 
These are calculated using the CSR-wake W CSR (z) and 
the bunch distribution X(z), where the average energy 
change per unit length (W CSR ) and the standard devia- 
tion aw (Wcsr) over the distribution are 



(W CSR ) = / Wosr^)A(^) dz, 



aw 



z + 



wl sn {z)\{z) dz-(w CSR y 



1/2 



(59) 



(60) 



The term aw is important because it contributes to the 
correlated energy spread in a bunch. 

To show how these quantities change as a bunch pro- 
gresses through a bend, Fig. [TT] plots (W C sr) and aw, 
normalized by Wq, versus different bunch center coordi- 
nates So in bend 1 using Eq. with D\ — lm and 
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Ki = K2 = 1/lOm. In the literature, the wake near the 
beginning of bend 1 is often calculated as if the prior drift 
length Di — > oo 0, so such calculations are plotted 
in dotted lines for comparison. From the difference be- 
tween the two approaches, one sees the effect of bend 2, 
where the CSR-wake at sq = is non-zero. In this ex- 
ample, they coincide after about 1.4m and 1.8m for the 
free space and shielded cases, respectively. 

In order for it to be plausible to ignore the vacuum 
chamber sidewalk, such a chamber must be wide enough 
to allow a straight path between the retarded bunch and 
the test particle. In this example, the vector from a 
source particle at z = —8000 a to the center of the bunch 
(z = 0) requires that the vacuum chamber half-width 
must be greater than approximately 3cm. 



for s' > B\. The total exit wake is then 



d£ c 



ds 



(s > Si) 



d£h 



+ 



d£ c 



ds 



ds 
dScsn 



D 



Bi 



(61) 



where the individual terms due to element elements Dq, 
Bi, Di, are written out in Appendix 1X1 For a magnet of 
length B\ = 3m, the exit wakes in the following drift Dq 
are shown in Fig. [T3] for bunch centers in the following 
3 meters between parallel plates and in free space. The 
average and standard deviation of the wakes through this 
region are shown in Fig. 1141 In the shielded situation, one 
sees that the bunch actually gains some energy in a short 
length following the bend, and that the total energy loss 
between parallel plates is negligible compared to the free 
space losses. Energy spread, however, is qualitatively the 
same in both cases. 



B. CSR in a Drift Between Bends 



The non-zero CSR-wake at the beginning of bend 1 
in Fig. [TT] is evidence that the wake in a drift region af- 
ter a bend also needs to be considered. This exit-wake 
in the region Dq following bend 1 is calculated using 
Eq. © with Eq. (J53J) and Eq. (JS7J around the center 
of a bunch at sq > B\. Because the bunch is moving 
in a straight line, the regularization procedure simply re- 
moves the need to integrate any s' > B\ for the real 
bunch. Therefore we can use Eq. for bend 1, drift 1, 
and earlier elements, and subtract the space charge terms 
for s' < B\. Image charges, however, still require terms 



V. BUNCH COMPRESSION 

Bunch compression or decompression can be achieved 
in a bending magnet if there is a correlation between en- 
ergy and longitudinal position of particles in the bunch. 
To exactly calculate CSR for this, however, requires at 
least a 2-dimensional model, because particles of differ- 
ent energies travel on different orbits. In the framework 
of the 1-dimensional model described by Eq. ^ , this ef- 
fect can be approximately modeled by allowing the bunch 
length to be time dependent, and neglecting variations in 
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-4-2024 -4-202 
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FIG. 13: CSR- wakes for various bunch centers so > Bi calculated using Eq. (|61[) , The left graph uses parallel plates separated 
by a distance H = 2cm, while the right graph is for free space (n = terms only in Eqs. (|A6|I - (|A7|) . and without Eq. l|A5|l V 
The bending radius kJ -1 = 10m, and the bunch has a Gaussian profile with a = 0.3mm and an energy of lGeV. 



the velocity /3 c. The density and current are then 

tt A na I r / s - s b - (3cb 

J(s,t) =Qf3c—X — 

a(t) V CT W 



(62) 



where A has unit norm and variance with respect to s, as 
in Eq (|3~7| . The time derivative of p(s,t) is 



d_ 

dt 



p(s,t) = -0c 



(7- 

V a 
a 2 



(63) 



with st = s— Sb~ (3 ct. Note that &/(/3 c) is on the order of 
a I B in a magnet of length B, and (s — Sb — (3 ct) is on the 
order of a for all relevant (s, t), and therefore the term in 
brackets is on the order of a/B <C 1 relative to the first 
term, and will be neglected. With such an approxima- 
tion, the CSR-wake in a bunch compression system can 
be modeled by simply making the substitutions 



A(s r ) 

A'(«r) 



1 



c(tret 
1 



-A 



X' 



t Iet = t - v/r 2 + (niJ) 2 /c) 



(64) 

(65) 
(66) 



in all of the previous formulas, with r = ||X(s) — X(s')|j, 
as in Eq. ([3]). This accounts for the real charges (n = 0) 
and image charges (n =/= 0) at the appropriate retarded 
times. 

Calculations for the average and standard deviation 
of the CSR-wake with a linearly compressing bunch 
through the length of bend 1 in free space are shown 
in Fig. 15(a) The above approximation is referred to 
as Method 1. Method 2 calculates the instantaneous 



CSR-wake of a compressing bunch at each point in the 
bend as if it always had its instantaneous length. Such 
a scheme is essentially what particle tracking codes (e.g. 
elegant [T3|, Bmad [HI]) use f° r CSR simulation. For 
reference, Method 3 calculates the CSR-wake for a non- 
compressing bunch that maintains the same length as the 
final compressed length in Methods 1 and 2. In this ex- 
ample, Method 2 overestimates the CSR effect compared 
to the more realistic Method 1 , and both exhibit a much 
smaller effect than Method 3. At the end of the magnet 
(sq = 3m), the CSR-wake, according to Method 1, has 
yet to reach its corresponding steady-state strength. 

Figure 15(b) shows these same calculations but be- 
tween parallel plates with H — 2cm. One sees that the 
energy loss in method 2 is similar to that in method 1, 
but the energy spread induced is overestimated. Free 
space and sh ielded calcul ations are repeated with D\ — > 
Figs. 15(c) -15(d) which when compared with 



Figs. [l5(a)}|T5(bJ one can see the effect of the previous 



bend B 2 - 



VI. COHERENT POWER SPECTRUM 

Some of the first CSR calculations are found in an 
originally unpublished report by Schwinger [9|. Here we 
use one of his methods to derive an exact expression for 
the coherent energy loss by a Gaussian beam, which is 
then used to verify our earlier calculations. Consider the 
power spectrum due to a single particle moving on a cir- 
cle with velocity /3 c, which is proportional to the absolute 
square of the Fourier transform electric field EW(fl,t), 
integrated over solid angle O, as in 



oo 

J dfl J dte lwt E (1) (r!,i) 



(67) 



For iV particles moving on this circle with positions s — 
s n + [3 ct, the total electric field can be written in terms 
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s (m) s (m) 

FIG. 14: The average energy loss and energy spread per unit length of the exit wakes in Fig. 1131 In this example, shielding by 
parallel plates drastically reduces the energy loss, but only marginally reduces the energy spread, when compared to free space 
calculations. 
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(a)Free Space Bunch Compression 
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(b)Parallcl Plates Bunch Compression 
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s (m) 

(c)Free Space Bunch Compression, D\ — > oo 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

s Q (m) 

(d)Parallel Plates Bunch Compression, D\ — » oo 



FIG. 15: The average and standard deviation of the CSR-wake in free space (Fig. 15(a) I and between parallel plates with 
H — 2cm (Fig. 15(b) l over a Gaussian bunch, compressing from a = 0.9mm to a = 0.6mm linearly through bend 2, and from 
a = 0.6m m to a = 0.3m m linearly through bend 1, using methods described in the text. Figures |15(a)p5(b)1 have Di = im, 
while Figs. 15(c) 15(d) have D\ — > oo. The lengths B\ = B2 = 3m, the bending radii are kJ -1 = k^" = 10m, and the energy is 
lGeV. 
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of the single particle's electric field (s n = 0), as in 



N 



EW(n,t) = ^EW(0,t-t n ), 



(68) 



n=l 



where the time deviations f„ = s n /(f3c). By changing 
variables, this means that the N particle power spectrum 
is simply 

2 



duj 



N 



dio 



(69) 



These phase factors can be separated into terms with 
m = n and m =/= n, 



/V 



du 



\m-l 



n=l 



N- 



iV 



A' 



(70) 



da; V /3c / z — ' V /3c 

so that the second term can be written as a correlation 
between different particles 



^exp »w- ^ rapl-iw— 



E ex P 



/3c 



x J ds \(s) exp w 
x y ds' A(s') exp 



N(N-1) 



pc 



(71) 



pc 

using the normalized particle distribution A(s) along the 
circle. The N particle power spectrum is then 



duj 



(to) ~ iV 




iV(iV - 1) 



JdsX( s )ex P (zg) 



dw 



(72) 



coherent 



The first term in Eq. (|72[) is the incoherent power spec- 
trum, while the second is the coherent power spectrum. 
The squared integral is called the form-factor. 

In free space, the well-known single particle power 
spectrum is 



(73) 



where u> c = §7 3 ck is the critical frequency [l0]. The 
function 5 is defined as 

CO 

S(0 = ^U / dxtf 5/3 (ar), (74) 
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FIG. 16: The power spectrum in Eq. (|72jl. per particle, using 
a Gaussian form factor with various values of the coherence 
parameter a c = 2 P , defined in Eq. (J77J). The lower frequencies 
are enhanced by a factor of N, and in this example N = 10 9 . 



in which K is a modified Bessel function. The integral 
J °° S{x) dx = 1, and the total power lost by a single 
particle is 



PW = -r c mc 3 /3 4 7 4 « 2 
3 



(75) 



For a Gaussian distribution with variance a 2 and n a <C 
1, the form-factor is, extending the integration limits to 
infinity, 



d.s 



WS 8 



27TCT 



exp 



~2 , ,2 



/3 2 c 2 



(76) 



: exp 



defining the coherence factor 



a c = - 7 ko 



(77) 



The total power spectrum per particle for an iV-particle 
Gaussian distribution with various values of a c is shown 
in Fig. [TBI One sees from the exponential that the lower 
frequencies, up to a cutoff frequency around u = /3 c/ a, 
are enhanced by a factor of N by the coherent part of 
Eq. (|72[) . The spectrum at higher frequencies agrees with 
the familiar single particle spectrum in Eq. (|T3[) . 

It turns out that Eq. (|72p can be integrated exactly 
for a Gaussian distribution. Explicitly, the total power 
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FIG. 17: The coherence function T c (a c ) of Eq. (J79J) is plotted 
in red. The green curve is the first term in the asymptotic 
expansion in Eq. (|80p [13], and the blue curve uses all three 
terms in Eq. ((8f])) . 



radiated by N particles is 

oo 

pM = jVpW f S(x) dx + N(N - 1)P (1) 



9y/3 



x e 



NP {1) + N(N - 1)P (1) 
9V3 



8tt 



K 5 /3{y) dy 



2 2 

~ a ^ x dx 



dx 



(78) 



dy 



= NPM + N(N - 1) P« T c ^ 7 3 k aj , 

in which the final integral yields the coherence function 
defined as 



T c (a c ) 



9 



32y / 7r a; 



i 



16a 2 



(79) 



The limit lim Q _, + T c (a) = 1, which is to say that an 
infinitely narrow bunch radiates as one charge. In prac- 
tical situations a c 3> 1, so an asymptotic expansion of T c 
gives the useful approximation 



T c (a c ) 



9T 



16 2 /3^Af 



1 



4/3 



9r 



32 • 2 2 / 3 ^ 



11 

16 \o c 

10/3 



(80) 



The first term in Eq. (JHOJ) is given in Ref. [II]. Fig- 
ure [T7| compares this first term to the exact expression in 
Eq. ([79]) and to all three terms in Eq. ([80]) . One sees an 
excellent approximation for a c > 50 using the first term 
and for a c > 1 using all three terms in Eq. ([50]) . Also, the 
average coherent energy lost per particle per unit length 
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FIG. 18: The relative difference \{b — a)\/b for a the average 
energy lost using Eq. (|32|l , and b being the result in Eq. (I81|l . 



is 



coh. 



r(f) 



VF + .. 



(82) 



using Wq defined in Eq. (|40p. The numerical coeffi- 
cient r(5/6)6~ 1 / 3 7r- 1 / 2 ~ 0.350. The same procedure 
in Eq. ([75]) and Eq. ([50]) can be carried out for a uni- 
form distribution of length AL with the same variance 
cr 2 , implying that AL = 2\/3a. The result yields the 
same form as Eq. ([82]) . except with the numerical coeffi- 
cient 2^ 4 / 3 ~ 0.397. This term was originally derived in 
Ref. [|. 

To verify that the CSR-wake does indeed represent the 
coherent energy lost, the relative difference of the average 
energy loss using the steady-state wake of a Gaussian 
bunch in Eq. ([32]) to the result in Eq. ([5T]) is plotted in 
Fig. [T5J One sees that the relative difference is at most 
1% in this practical parameter range, and that occurs 
with relatively long bunches. We speculate that this error 
is caused by the regularization procedure that subtracts 
the space charge term from the longitudinal electric field. 

The relevance of the coherence function depends on 
the number of particles N — 1 ~ N. The coherent power 
radiated equals the incoherent power radiated when iV • 
T c (a c ) = 1, illustrated in Fig. [19] Using Eq. (J80j, the 
coherent power dominates the total power when 



< 



AT3/4 



(83) 
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FIG. 19: The dividing line where the coherent power equals 
the incoherent power, i.e. the total power is twice the inco- 
herent power. Below this line, the coherent power dominates 
the total power. 



VII. CONCLUSION 

The wake-field due to CSR of a 1-dimensional bunch 
traveling on a curve without small angle or high energy 
approximations has been derived using Jefimenko's forms 
of Maxwell's equations. This exact solution allowed us 
to quantify the accuracy of the approximations of the 
steady-state CSR-wake in a bend given in Ref. 4] and 
Ref. [7| showing that the former is inaccurate at low en- 
ergies and long bunch lengths, and that the latter is much 
more accurate down to low energies. All approximations 
tend to overestimate the CSR-wake. For planar orbits 
the equations are extended to include shielding by per- 
fectly conducting parallel plates using the image charge 
method. 

The formulas have been applied to the geometry of 
a bend preceded by a drift, preceded by another bend, 



and show that the CSR-wake well inside the downstream 
bend is influenced by the upstream bend for the param- 
eters used. In fact, a bunch near the entrance of a bend 
is influenced by the CSR-wake due to the previous bend 
much more than by that due to the previous drift. Shield- 
ing by parallel plates reduces the energy loss rate signifi- 
cantly, but the effect on reducing energy spread increase 
is far less dramatic, in both the drift and bend regions. 

Bunch compression has been added to this model by 
allowing the bunch length to be time dependent, so that 
the retarded charge density seen by a test particle is ap- 
propriately taken into account. This method has been 
compared to simple methods used by particle simulation 
codes Bmad and elegant, and it is shown that these tend 
to overestimate the effect [IJ B3 ■ 

Additionally, an exact expression for the coherent 
power lost by a 1-dimcnsional Gaussian bunch moving in 
a circle has been derived by integrating the power spec- 
trum, following the method of Schwinger [9( • When com- 
pared to the energy loss rate by the CSR-wake, the two 
show slight deviations. This could be due to the regular- 
ization procedure for the 1-dimcnsional CSR-wake that 
subtracts off the space charge term. 



APPENDIX A: FORMULAS FOR MULTIPLE 
BENDS AND DRIFTS 

For ease of reading, the individual terms terms in 
Eq. ([58)1 and Eq. ([rjT)l have been deferred to here. 
They are calculated by applying Eq. ([3]), regularized by 
Eq. (fTTJ|) . including image charges as in Eq. ([4"5)l . to the 
geometry in Eq. ([53f . 

In Eq. (|58|) . the first term d£ csr /ds \b ± is the sum of 
Eq. ([29| and Eq. ([47)) with k — ► K\ and 9 — ► K% s, explic- 
itly 



d£ cs 
ds 



■<s) 



-Nr c mc 



/ 



da 



j3 2 cos (a) — 1 1 sgn(a) — j3 cos (a/2) 
2|sin(a/2) | + a - 2/3 1 sin (a/2) | 



A'(« a ) 



Kl X(s a ) 


oo 

at <• 

+ / dA 


1 X'(z- 


2| sin (a/2) | 


7 2 A 




A a 




OO 


-«1 X(s a . n ) 


+ 




>ot,n 





dA 



1 A'(z + A) 



(Al) 



da 



P 2 cos (a) — 1 , 



A (s a ,«) 
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with the definitions 



a a = Ki(s- Bi), 
at = Ki s, 



Aa 



r Q , n = — 2 cos a + (n K\H) 2 , 
s a = s — so (a — /3\/2 — 2 cos a) , 

Kl 

= s - s (a — r an ) . 

Kl 

Some trigonometric functions have been simplified, and the space charge integrals have changed variables to A = 
(a — 0\J 2 — 2 cos a)j k\. These terms account for the regularized CSR-wake and image charges in bend 1. The next 
terms are 



D, 



ds 



= Nr r mc 



f dLjr(2-£„ >0 )(-iy 

n n=0 



R 



, cos [k\ s) T L 



Rh,n 



R 2 

L, 



X'(s L , n ) 



Rl.u = — \/2 — 2 cos («i s) + 2k\L sin (k\ s) + (k\L) 2 + (kitiH) 2 



Tl = L cos (ki s) H sin (k\ s) 

SL,n = —L — s + (3R L ^ n 



(A3) 



ds 



= Nr r m& 



/ dL^(2-<5„, )(-l) 

n n=Q 



R 



A(s L ,„) + 



L,n 



, cos (ki s + k 2 L) Tl 



Rh,n 



-0 



R l,n 



X'(s L , n ) 



Rh,: 



COS {K\ s) — 1 1 — COS («2 £) 



«2 



£>i 



sin (ki s) sin (K2 L) 



Kl 



K2 



(nH) 2 



(A4) 



Tl = D\ cos (ki s) H sin (ki s) H sin (ki s + k 2 L) 

KlK 2 k 2 

Si )Tl = -L - Di - s + R L , n - 

Note that the lower limit of the sums have been set to n — to account for the real charges as well as image charges, 
necessitating the use of Kronecker's delta. Alternatively, if only free space terms are desired, the above formulas can 
be used with the n = term only. The dummy variable s' has been rescaled to L which integrates backwards over 
the length of the appropriate element. The terms Rl,u, Tl, and sl,u are redefined after each equation in order to 
keep the naming sane. 

Similarly, the wake at s > B\ after bend, as in Eq. (|6ip . contains the terms 



d£ c 



ds 



-Nr c mc 




, 1 m A(s' - s + 0(s - s')) + 0^—^X'(s' - so + 0(s - s')) 
(s — s 1 Y s — s 



XjB.-so + ^is-B^ + jnH) 2 ) 
h ^/(s-B^ + inH) 2 



(A5) 



OO 

OO n 

■X)2(-l) n / 

n=l { 



dL 



X'(L + Ih- sp + /V(s -B Y - L) 2 + (nH) 2 ) 
7 V(s-^i-^) 2 + (nH) 2 
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d£ c 



ds 



Bi , 

- Nr c mc 2 £(2 - «J n ,o)(-l) n / dL \ * A(s L ,„) 

n=0 { { K L,n 



Rl,t, 



2 cos(kiL) _ Tl_ 

P r> P 7->2 

L.n 



Rl,u R 



A'(SL,«) 



12 2cob( Ki L) + ^-BQainfoS) + (g _ ^ + (nff)2 



K 2 



Kl 



T L = s - Bi H sin(/ciL) 

SL,n = -i + Si - S + P RL.n, 



(A6) 



d£ c 



ds 



= Nr c mc 2 J2(2-5n.o)(-l) n / dL J * A(«£,„) + 



2 COS_(Kigi) 
Rl,u 



R l.n 



\\s L ,n) 



Rl.u = \ \ L+{s-B 1 ) cos (kiBi) + 



sin(Ki£?i)\ /cos (/ti-Bi) — 1 



Kl 



+ 



Kl 



(s-Bi)sin(KiBi) +(nfl-r (A7) 



Tl = s — Bi + L cos (ki-Bi) H sin (ki-BI) 

Kl 

= —L — so + 0y/Rl + {nHf. 
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